Disassortativity of random critical branching trees 
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Random critical branching trees (CBTs) are generated by the multiplicative branching process, 
where the branching number is determined stochastically, independent of the degree of their ancestor. 
Here we show analytically that despite this stochastic independence, there exists the degree-degree 
correlation (DDC) in the CBT and it is disassortative. Moreover, the skeletons of fractal networks, 
the maximum spanning trees formed by the edge betweenness centrality, behave similarly to the 
CBT in the DDC. This analytic solution and observation support the argument that the fractal 
scaling in complex networks originates from the disassortativity in the DDC. 
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Recently, it was discovered [T| that many complex net- 
works in real world are fractals, satisfying the fractal 
scaling: The number of boxes Nb(£) needed to cover 
an object scales in a power-law manner with respect 
to the box size I, i.e., Nb{1) ~ t~ dB , where ds is the 
fractal dimension. Examples are the World-Wide Web 
(WWW) 0, the protein interaction network (PIN) of 
budding yeast [|| and the metabolic networks In 
contrast, the Internet || and many artificial model net- 
works such as the Barabasi- Albert (BA) model @ and 
the static model Q, etc, are not fractals. It was argued 
that the fractal scaling originates from the disassortative 
correlation between two neighboring degrees [|[ or the 
repulsion between hubs @. 

The origin of the fractal scaling has been understood 
from another perspective (lolllH: A network is composed 
of the skeleton, which is a special type of spanning tree 
formed by edges with the highest betweenness centrali- 
ties or loads, and the remaining edges in the network that 
contribute to loop formation. For fractal networks, it was 
shown that the skeletons exhibit fractal scaling similar to 
that of the original network. The number of boxes needed 
to cover the original network is almost the same as that 
needed to cover the skeleton. Moreover, when a skeleton 
is considered as a tree generated in a branching process 
starting from an arbitrary selected root vertex, the mean 
branching number, the average number of offsprings, ex- 
hibits a plateau, albeit fluctuating, independent of the 
distance from the root. The value is close to 1, and 
the skeleton was regarded as the critical branching tree 
(CBT), which is known to be a fractal [121] . Thus, the 
fractal scaling in the original network originates from the 
presence of the fractal skeleton underneath the original 
network. 

The CBT is generated by the multiplicative branching 
process. To generate a scale-free tree, n (> 0) offsprings 
are created at each branching step with the probability 
b n) which is given as follows: b n = n 1 /((j— 1) for n > 1 
and bo = 1 — Xm=i where Ci x ) is the Riemann zeta 
function. Then the obtained branching tree is a scale-free 
network with degree exponent 7. Since branching event is 
stochastically independent, one may think that the CBT 



is random in the degree-degree correlation (DDC); how- 
ever, here we show that the DDC is disassortative. We 
also show that the skeletons of the fractal networks also 
exhibit the similar mixing pattern. Therefore, the origin 
of the disassortativity of the fractal networks is rooted 
from the CBT nature of the skeleton. 

Here, we calculate the two point correlation function 
P(k, k') for the CBT. P(k, k') of an undirected network 
is defined as the fraction of links with degrees k and k' on 
both ends. Even though the network under consideration 
is undirected, for further discussion, we make it directed 
by assigning arrows to each link in an arbitrary manner. 
Then we count the number of links with degree k on 
the arrow's source side and k' on its sink side and call it 
N(k — > k'). Next reverse all arrows of the links and count 
the same and call it M(k — » k'). Each link contribute 
once in N(k -> k') and M[k -> k'). Then 



P(k,k') 



N{k -> k') + M{k -> k') 
2L : 



(1) 



where 2L = (k)N is twice of the link number. Note 
that M(k -> k') = N(k' -> k). This way, a (3-1) link 
contributes to the element P(l,3) once and P(3, 1) once 
while a (2-2) link contribute to P(2, 2) twice. Since the 
sum is normalized by 2L, we have the general relation 



J2P(k,k') 



kP d (k) 
(k) 



(2) 



with Pd(k) the degree distribution of the net- 
work. For uncorrelated networks only, P(k, k') — 
kP d (k)k'P d (k')/(k} 2 0. 

Using above procedure, P(k, k') for the CBT with 
(k) = 2 is obtained as follows: First, we consider a large 
enough CBT so that we may neglect the boundary effect. 
Assign arrows in the natural way following the branch- 
ing direction. Then N(k — > k') is the number of degree 
k nodes (NP d (k)) times the number of offsprings (k — 1) 
times the probability that those offsprings has degree 
k'. So, A^(fc -» k') = NP d (k)(k - l)P d (fc'). Similarly, 
M(k -» k') = NP d {k'){k' - l)Pd(fe). So, we find 



1 



P(k,k') = -(k + k'-2)P d (k)P d (k'), 



(3) 
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FIG. 1: (Color online) Plot of (fe nn )(fc) as a function of degree 
k for the CBT (o) and the skeletons of the fractal networks, 
the WWW (□) and the yeast protein interaction network (A). 
Solid lines represent the formula, Eq. ((5j). Inset: To visual- 
ize the 1/k dependence, y(k) — ((k nn )(k) — a)/b is plotted. 
Eq. (0 predicts y(k) = 1/k, represented by the solid line. 
The numerical data fit this form reasonably. 



which is different from the uncorrelated ones, counterin- 
tuitively. 

The DDC manifests in the mean degree (k nn )(k) of the 
nearest neighbors of a node with degree k. It is related 
to P(k, k') as 



{k nn )(k) = Y / 



k'P(k,k') 
kP d {k)/{k)' 



and is independent of k for uncorrelated networks, 
ging the formula §5§ into Eq. we obtain that 



<fcnn>(fc) = 



(fc) 2 (fc)((fc 2 )-2(fc)) 
2 2k 



(4) 



Plug- 



(5) 



(k) = 2 for the CBT. Eq. (J5j) may be rewritten in 
the form, (k nn )(k) = a + b/k, where a = {k) 2 /2 and 
b = [(k)((k 2 ) - 2{k)))/2. Thus, (k nn )(k) is inversely pro- 
portional to degree k for k < b/a and thus the CBT is 
disassortative. 

We check this disassortative behavior numerically for 
the CBT with 7 = 2.5 averaged over 10 samples with size 
N = 10 6 and show it in Fig. [T] Indeed, the numerical 
data fit well to the analytic result Eq. ([5]), represented by 
the solid line in Fig. 1 and its inset. We note here that for 
7 < 3, b/a is large and scales with N as ~ A^ 3-7 ^ 7 " 1 ), 
which we confirm numerically. 

Next we consider the skeletons of fractal scale- free net- 
works. As discussed above, we argued they could be ap- 
proximated as CBTs. To corroborate it, we also show in 
Fig. 1 (fc„„)(fc)'s of the skeletons of the WWW and the 
PIN and compare them with Eq. (5) shown as solid lines 
where the measured values of (k 2 ) are used together with 
(k) = 2. We find the agreements quite good. Thus, these 
skeletons can be regarded as having the same DDC as 
that of the CBT. On the other hand, the skeletons of non- 
fractal networks show different behaviors. Although the 



Internet at the autonomous system (AS) level exhibits a 
disassortative mixing pattern, the decaying behavior of 
(k nn )(k) for its skeleton is different from that of the CBT 
as shown in Fig. [5] It decays as (k nn )(k) ~ fc -0 7 . For 
an artificial model, e.g., the static model with degree ex- 
ponent 7 w 2.4 < 3, (k nn )(k) of its skeleton decays as 
~ fc~ - 8 as shown in Fig. [51 different from ~ fc _1 for the 
CBT. For the BA model with degree exponent 7 = 3, 
{k nn )(k) is random; however, for its skeleton, it is weakly 
disassortative for intermediate range of k. 




FIG. 2: (Color online) Plot of {k n n)(k) as a function of degree 
k for the skeletons of the non-fractal networks, the Internet at 
the AS level (o), the static model with the degree exponent 
7 = 2.4 (A), and the BA model with the degree exponent 
7 = 3.0 (□). Solid (dashed) line is a guideline with slope 
-0.7 (-0.8). 

In summary, we have shown analytically and numer- 
ically that the CBT is disassortative in the DDC. This 
is induced topologically through the branching process. 
Its origin is similar to what was proposed [l5[ and shown 
later [Ty] for the Internet that the disassortative mixing 
pattern is caused by the topological restriction that no 
pair of nodes is allowed to have multiple connections in 
the ensemble of graphs with given (or expected) degree 
sequences. In the CBT, such restriction arises among the 
offsprings of a same ancestor that cannot be connected 
each other. Here, we also showed that the skeletons of 
fractal complex networks such as the WWW and the PIN 
exhibit the same pattern in the DDC as found in the 
CBT. This is yet another evidence that the skeletons of 
the fractal networks can be regarded as the CBTs besides 
the mean branching ratio b eing close to 1 independent of 
the distance from the root jldl [TlT | . On the other hand, 
the skeletons of the non-fractal networks show different 
patterns in the DDC, even though they are disassorta- 
tive. The explicit formula ([3]) derived here can be used 
to study various dynamic problems on fractal networks 
such as the epidemic problem and so on. 
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